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We determine the full post-Newtonian limit of theories of gravity that extend general relativity 
by replacing the Ricci scalar, R, in the generating Lagrangian by some analytic function, f{R)- We 
restrict ourselves to theories that admit Minkowski space as a suitable background, and perform a 
perturbative expansion in the manner prescribed by the parameterised post-Newtonian formalism. 
Extra potentials are found to be present that are not accounted for in the usual treatment, and a 
discussion is provided on how they may be used to observationally distinguished these theories from 
general relativity at the post-Newtonian level. 
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■ I. INTRODUCTION 

K-5 There exists an extensive literature on relativistic theories that generalise Einstein's theory of general relativity 

J — . (GR), and that reduce to GR in the appropriate limits. A particularly appealing class of these generalisations are the 
fourth-order theories. These are theories derived from a Lagrangian density that is a scalar function of contractions 
of the Riemann tensor only. Considerations of fourth-order theories of gravity have a long history, having been first 
considered by Eddington in as early as the 1920's [Tl|. One frequently considered generalisation is to replace the Ricci 
curvature scalar, R, in the Einstein-Hilbert action with some analytic function, f{R) (see e.g. 0,0, 13: Q)- It is the 
post-Newtonian limit of such theories that we will be interested in here. 

There are a variety of reasons why one may wish to consider these generalised theories. Strong motivation comes 
the from their renormalization properties in the presence of matter fields Other motivation can be found from 
cosmological considerations where it has been found that generalising the Einstein-Hilbert action can be of use for 
. better understanding the late-time acceleration of the universe 0, 111 , early universe inflation 0, [H, [ll| and the nature 
' of the initial singularity [I2j IH, Hi [III ■ Whatever the motivation for considering generalised fourth-order theories, 
it is essential that we sufficiently understand their weak- field limit, and that they conform with the ever increasing 
' body of observational data. 

The usual frame- work for considering the weak-field effects of modified theories is the parameterised post-Newtonian 
(PPN) formalism. The first steps towards this formalism were again made by Eddington T| , together with Robertson 
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^0 ', tlfil and Schiff [17[, who treated the planets as test bodies moving in the gravitational field of the sun. In their 



formalism they introduced the test metric 



,9 (2m 2Bm^\ / 27771 \ , ^ 

' where dx^ is the three dimensional Euclidean line-element, r is a radial coordinate and 777 is the mass of the central 
gravitating object. The parameters /? and 7 are to be determined by experiment. By constraining them it is possible 
to verify, or potentially disprove, GR and its alternatives. For example, GR predicts [3 — ^ — 1, whilst Brans-Dicke 
theory [ISj predicts /3 = 1 and 7 = (1 -I- w)/(2 -I- uj), where uj is the Brans-Dicke coupling constant. Using radio 
communications from the Cassini spacecraft the constraint 7 = 1 -f (2.1 ± 2.3) x 10"^ has been obtained [l^. This is 
evidently in good agreement with GR, and can be used to constrain the Brans-Dicke parameter to be a; > 40 000, to 
2cr. 

Our goal here is to calculate the post-Newtonian limit of fourth-order theories, so that they can be the subject of 
observation in a similar manner. To achieve this we will consider them in the context of the full PPN formalism. 
This formalism is a generalisation of the Eddington, Robertson, Schiff parameterisation outlined above, and allows 
for more general configurations of matter fields than a single point-like gravitating source. The PPN approach was 
developed primarily by Nordvedt and Will ^2^, and is explained in detail in ^T^. We will give a brief explanation of 
the PPN approach below, in so far as will be required for coherence of this work. For a more complete exposition the 
reader is referred to I21I. 
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Absent in previous studies is a rigorous exemplification of how fourth-order theories fit into the PPN formahsm. 
Here we will remedy this by extending the PPN formalism to include fourth-order /(i?) theories. In doing so we will 
attempt to maintain, to the highest degree possible, the principles and spirit of the PPN formalism, as expounded in 
[2l| . Previous attempts have been made in this direction by Capozziello and Troisi [l^l and Olmo J23\ . These authors 
attempt to derive the post- Newtonian limit of f{R) theories by appealing to their equivalence with scalar-tensor 
theories [ISl- Here we work directly with the fourth-order theory, and by direct integration of the field equations 
find results that explicity state their post-Newtonian limit. In doing so we restrict our attention to the subset of 
theories that admit Minkowski space as a suitable background, and find that it is necessary to introduce a number 
of new post-Newtonian potentials. Theories with other backgrounds will be investigated elsewhere. We note that the 
post-Newtonian limit of Gauss-Bonnet gravity has been found by Sotiriou and Barausse [1^ . 

In section [H] we introduce theories that are derivable from a Lagrangian density of the form C = f{R). The field 
equations are derived, and a discussion is given of the condition that Minkowski space be an appropriate background. 
Section [Till gives a brief introduction to the PPN formalism, in as much as is required for the self-consistency of this 
article. In sections ITVl and fVl we find the Newtonian and post-Newtonian limits, respectively. These calculations are 
performed in the presence of a perfect fluid, and result in many new potentials that are not usually present. In section 
IVII the results found in the previous two sections are transformed into the standard post-Newtonian gauge, in which 
the spatial part of the metric is diagonal and derivatives of quantities associated with the matter fields are removed. 
Section [VIII gives a discussion of the results obtained, and in particular gives the relevant post-Newtonian limit if the 
usual Newtonian potential is to dominate at the Newtonian level of approximation. In section rVIIII we conclude, and 
indicate the principle ways in which the f{R) theories we are considering may be observationally distinguishable from 
GR at the post-Newtonian level. The appendices give some details of the more lengthy calculations. 



II. FOURTH-ORDER THEORIES 



The Lagrangian density for the theories we will be considering is 

C = f{R). (1) 
The action associated with this density is then given by 

J V^C + S^, (2) 

where Sm denotes the action associated with the matter fields. Extremizing this action with respect to the metric 
results in the field equations 

f'Rf^i^ - ]^f9piu + f'-^^{g^ugap - 9nag,,p) = SttT^^, (3) 

where primes denote differentiation with respect to R and T)„y is the energy-momentum tensor, defined in terms of Sm 
and gfj^iy in the usual way. Throughout we us Greek letters to run over all space-time indices and choose units so that 
c = G = 1. The Lagrangian formulation of these theories guarantees the covariant conservation of energy-momentum. 

In order to define a perturbative expansion we must first decide the appropriate background to expand about. In 
the usual PPN treatment this background is taken to be Minkowski space, and the metric is then expanded as 

gpi- = Tlf^f + (4) 

where h^,^ <^ 1. Such an expansion is well motivated in GR where it is known that in the absence of any matter 
fields Birkhoff 's theorem ensures staticity of spherically symmetric space-times, and that these space-times will be 
asymptotically flat. Adding small amounts of matter to such a background is then well modelled by the perturbative 
expansion given by ([4]), and one need only be concerned with the effects of matching this region to a suitable 
cosmological solution at large distances. 

More care is required with fourth-order theories, where staticity and asymptotic flatness cannot be so easily assumed. 
BirkhofF's theorem is not valid in these theories, and so staticity, if it is required, must be imposed as an extra condition 
on spherically symmetric vacua. However, the imposition of this extra symmetry does not, in general, result in 
Minkowski space as a suitable background. It was shown explicitly in [2^ that if any single power of the gravitational 
Lagrangian dominates, other than the Einstein-Hilbert one, then spherically symmetric vacuum space-times are not 
asymptotically attracted to Minkowski form. Instead, the line-element asymptotically approaches 
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as r — > oo, where n = only if the Einstein-Hilbert term dominates. Such asymptotic behaviour can be readily shown 
to be incompatible with observations, unless n is very small ,26J . 

An alternative prescription to the imposition of a static background was investigated in [23| ■ Here the method of 
imposing time independence, as outlined above, was contrasted with the method of assuming asymptotic homogeneity 
and isotropy of the background vacuum. The assumption of homogeneity and isotropy removes the need for r 
de pen dence in the asymptotic form of the metric, but at the expense of introducing time dependence. It was shown 
in [27| , using exact solutions as well as perturbative expansions, that the choice of symmetries for the background has 
a demonstrable effect on the weak-field expansions that are performed within them. 

In short, one must make a choice of background to expand about, and this choice can have important consequences 
for the expansion itself. If any single power of the gravitational Lagrangian other than the Einstein-Hilbert one 
dominates at asymptotically large distances, then Minkowski space may not be an appropriate, stable choice of 
background. One is then forced to recant either the time independence or the homogeneity of the background metric, 
with non-trivial consequences. 

Here we will avoid these difficulties by considering only theories in which the Einstein-Hilbert term dominates in the 
low curvature regime. We are then justified in performing a perturbative expansion about Minkowski space, which 
allows a more direct comparison with the usual PPN approach. The post-Newtonian limit of other f{R) theories 
will be dealt with in a future study. Fourth-order theories that admit a Minkowski background, and are an analytic 
fmictions of R, can then be written as 

oo 

/(i?) = 5]c,i?^ (5) 

2=1 

where the Ci are a set of real, positive valued constants. The post-Newtonian analysis that is to follow will show 
how the Ci are manifest in the weak- field limit, and hence how they can be potentially observed with gravitational 
experiments and observations. 



III. THE PPN APPROACH 



This section is a recapitulation of the PPN formalism, as propounded in [21|, and as is necessary for coherence of 
this article. The PPN formalism is a perturbative treatment of weak-field gravity. Such an expansion requires a small 
parameter to expand in. An "order of smallness" is therefore defined by 

f/o.t;2_^^n - 0(2) 
P 

where U is the Newtonian potential, v is the velocity a fluid element, p is the pressure of the fluid, p is its rest-mass 
density and H is the ratio of energy density to rest-mass density. Time derivatives are also taken to have an order of 
smallness associated with them, relative to spatial derivatives: 

\d/dx\ ^ ' 

(Recall that we have chosen to set c = 1). The PPN formalism now proceeds as an expansion in this order of smallness. 

The equations of motion show that for time-like particles propagating along geodesies the level of approximation 
required to recover the Newtonian limit is goo to 0(2), with no other knowledge of the metric components beyond 
the background level being necessary. The post-Newtonian limit for time-like particles requires a knowledge of 



500 to 


0(4) 


goi to 


0(3) 


gij to 


0(2). 



Latin letters are used to denote spatial indices. To obtain the Newtonian limit of null particles we only need to know 
the metric to background order: Light follows straight lines, to Newtonian accuracy. The post-Newtonian limit of 
null particles requires a knowledge of 500 and gij both to 0(2). 

Using the expansion ^ we can now calculate the Ricci and energy-momentum tensors to the appropriate orders. 
However, before doing so it is worth recognising that we have four gauge freedoms, associated with four coordinate 
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choices. Specifying the four gauge conditions 



hio,i = ^hiifi + 0{5) (6) 
^yj = ^f^jj,^ - ^f^ooj + 0(4) (7) 
allows the components of the Ricci tensor to be written 

Roo = -^V^/ioo - ^|V/iooP + ^hjkhoo,jk + 0(6) (8) 

i?o. = -^V2/io,-iv,o + 0(5) (9) 

i?y =-iv2/iy +0(4) (10) 

where — d^di is the Laplacian on three dimensional Euclidean space. We still have the freedom to make gauge 
transformations of the form 

and we will use this freedom in the following analysis to transform to a "standard post-Newtonian gauge" in which 
the spatial part of the metric is diagonal, and terms containing time derivatives are removed. The components of the 
stress-energy tensor, to the relevant order, are 

Too -p(l + n + i;2-/ioo) (11) 
Toi = -pvi (12) 
Ttj = pvivj +p6ij. (13) 

We can now substitute these expressions for i?^^ and Tfj_i, into the field equations ([3]), together with (O, and solve 
the equations order by order in perturbations. Transforming to an appropriate gauge will then yield the PPN limit 
of these fourth-order theories. 

IV. THE NEWTONIAN LIMIT 

The Newtonian limit of these theories will now be investigated. This limit has been found before for a point-like 
mass at the origin, originally by ^28] , and again several times since. Here we find the general solution for a space-time 
containing a perfect fluid. 

Beginning with the trace of the field equations ([3]), we have to 0(2) 

V^i?(^)--^^(^)=-i^p (14) 

6C2 3C2 

where R^^^ denotes the Ricci scalar to 0(2). This is an inhomogeneous Hclmholtz equation which has the solution 

3c2i |x-x'| ^ ' 

We have ignored here the other possible root in the exponential, which is an equally valid solution of the Helmholtz 
equation, but does not give an appropriate limit at asymptotically large distances. It can be seen that (jlSp is an 
exponentially decaying Yukawa potential if ci and C2 have the same signs, and is a damped oscillatory function if they 
have opposite signs. 

The — and trace field equations, given by ([3]), can now be written to 0(2) as 



Qci/.(^„)-Kicr/.r+2c2i?(^))=-8.p. 

and 

647r 



V 
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where the 0(2) parts of ([8]) and (llip have been used. Solving these two inhomogeneous Poisson equations simultane- 
ously we find 

C = -(C^ + C2i?('^) (16) 

where 

is the Newtonian potential. From pB]) it can be seen that the Ricci scalar itself acts as a Newtonian level potential. 



V. THE POST-NEWTONIAN LIMIT 



In this section we investigate the Post-Newtonian limit of the these theories. Expressions for hij to 0{2), hoi to 
0(3) and Hqq to 0(4) are found. These expressions will not be in the standard post-Newtonian gauge, discussed 
above. We perform a transformation into this gauge in the subsequent section. 



A. The hij terms 

We begin by evaluating the terms hij to 0{2). These quantities, together with (fT6|) above, are sufficient to determine 
the post-Newtonian limit of null geodesies. 

The i — j field equation ([3]) can now be written 

y z Ci Ci J 

where we have made use of the expressions above for i?^ and Ty, (fTO|) and ([T3| . the trace equation l|14p and the 
definition of U, ([TT]). This equation can be integrated to give 

= ^ - C2%i?(2) - 12|i?(5 + Afu,,>j . (18) 

Equation (jlSp is not diagonal, and so is not in the standard post-Newtonian gauge. In the next section will remove 
the off-diagonal components with the appropriate transformation. 



B. The hoi terms 



As discussed above, the first non-zero contribution to the hoi terms is at 0(3). The — i field equation ^ can now 
be written 

where we have used the expressions for Roi and Toi, © and (fT^ . the expression for h^o to 0(2), ([TB|) . the trace 
equation, p4|) . and the two new potentials Vi and Wi, that are defined as 



W. . / P(xO(v(x-) -(x-xO)(.-..0.^3,, (20) 



as in the usual PPN treatment, so that V (W^i — Vi) — 2Um- Use has also been made of the conservation equation 

|^ + V.(pv)=0. (21) 



Integrating the field equation above we now have 

~-^V.- ^W. + lO^C/,0, - 3o4 ^^'^ 
2ci 2ci cf cf 

which, again, will be subject to a gauge transformation in the next section. 



^ = - — ^» - — + lO^t^.O' - (22) 
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C. The hoo term to 0(4) 

At this order of perturbation the equations become more unsightly, and so we choose to relegate the majority of 
them to appendices, stating here only the results. Hence, in the gauge specified by (O and ([7]), the /iqo term to 0(4) 
is 



244lVi?|- 



7 



ISttci 
2 



V(C/i?) 



3C2 

47rci 



+ -V(pn) 

Cl 
8C2 



Cl 



-V(P) - -L (£1 - iL) X(R=) + ^X{UR) 

Cl 47r \ Cl 2C2 / DTTCl 



8C2 

3c2 



8c 



3c? 



Cl 



3ci 



00- 



The derivation of this result can be found in appendix [XI where it is given there as equations (|A3p and 
new potentials V, X and x £^re defined as 



V(Q) 



Q' 



-dx 



/3 



|x-x'| 



Q'e 



X' 



(23) 
The 

(24) 

(25) 
(26) 



Primes here label quantities that are functions of x'. It should be noted that this definition of V is degenerate with 
some of the usual PPN parameters: For example, V{pv'^) is identical to the potential <I>2 of [2l|. We use this definition 
of V as it is convenient for expressing the new potentials. 

Again, this result is not in the standard post-Newtonian gauge. In the following section we will transform it so that 
the terms proportional to i?,oo, t^,00: |VC/p, |Vi?p, VU ■ Vi? and x.oo are eliminated. 

VI. GAUGE TRANSFORMING 

In the preceding section we used the gauge specified by conditions ^ and ([7]). This has been a convenient choice, 
and has allowed integration of the field equations to post-Newtonian accuracy. However, it is desirable to transform 
the results found above to a gauge in which the spatial part of the metric is diagonal, and in which the metric takes it 
simplest form. By making the coordinate transformation x^ ^ x^^ + ^'^ the metric is transformed in such a way that 



^ ^l^':^ ^(C )• 



Then by making the choices 



Co 



6^C/.o 



18— T-i? 



6c?^'° 



the metric perturbations transform as 



,(2) 



24^i?,,, 



,(3) 



''00 



,(3) 




C2 I 
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(2) 

whilst /iqq is unchanged. These transformations are exactly what is required to diagonalize the spatial part of the 
metric, and to remove unwanted terms from the other metric components. The final form of the perturbed metric 
can now be written to the required order as 



1 (t/ + c,R) -^U' + 24i?^ - ^C/i? - T^nUR) + ^V{R') + 
ci Ci 3ci ISttci 4ttci 9c^ 



- ^V{pR) ~ f?V(Vp • Vt/) + ^V(Vp • Vi?) + -V(pn) + -V{pv') + -V{p) 

- ^X(Vp • Vi?) - -X{p) + — X(pn) (27) 

go. = - — - — + ^^^'''^ - ^ - (28) 
2ci 2ci 6ci 6ci 6V6ciC2 

5y = (^1 + ^ ([/ - C2i?)^ % (29) 
where we have introduced the new potentials Yi and Z^, which are defined as 



_ f pV • (x-x')(a;-x')j^-^/^|x-x I 3 , 



Y^^j - J— -' e V - 'd^,' (30) 



and where use has again been made of the conservation equation l|2ip . The reader will notice in equation (j27p 
potentials that are functions of gradients of p, U and 7?, such as V(V/o • VC/). This type of term is not of the usual 
PPN form, where the metric contains functionals of rest mass, energy, pressure and velocity, but not their gradients 
PH . In appendix [Bl we re-express the offending terms in a more proper PPN form, where gradients are absent. 



VII. DISCUSSION 



We have found the PPN limit of analytic /(i?) theories of gravity that allow an asymptotically Minkowski back- 
ground. The weak-field metric for these theories, in the presence of a perfect fluid and in the standard post-Newtonian 
gauge, is given by equations ((27)) . ((28)) and (|29|) . We shall now proceed to investigate their form. 

Firstly, we will consider the limit where the higher order contributions to the action are vanishing, so that C3 
C2 0. In such a limit it can be seen from the trace equation, ([U]), that R %T^p/ci, and that the metric specified 
by 1271), (UHl) and ^ then reduces to 

<7oo ^ - 1 + 2C/ - 2C/2 + 4V(pC/) + 2V{pIl) + AV{pv^) + 6V(p) 
7 1 

g^j ^ (1 + 2U) S^j 

where the two terms V{pU) and V{RU) in (j27p have contributed to V{pU) in the expression above. Here we have set 
ci = 1. This metric is the PPN hmit of GR (see [2l| for details). 

In order to be considered viable for non-zero C2 and C3 it is necessary for these theories to reduce to Newtonian 
gravity in the appropriate limit. From the 0(2) term of the goo component of the metric, (j27p . it can be seen that 
this can occur iff either 

{i) 3c2R - [/ or (ii) 3c2R ~ 

on observable length scales. Condition (i) is met if 

|Vci/6c2|L< 1 
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for the largest length-scales on which Newtonian gravity has been observed, L. Alternatively, condition (ii) can be 
met if i? is a decaying exponential with 



VcJGc^I > 1, 

where I is the smallest length scale on which Newtonian gravity has been observed. In both cases we can now find 
simple expressions for the post-Newtonian limit, which will allow us to relate our results to the relevant observations. 

Let us first consider case (i). This case can be easily dismissed by considering only the Newtonian limit of goo and 
the post-Newtonian limit of gij, equations (I27p and (j29|l above. If ■\/ci/6c2|x — x'| ^ 1 then 3c2R ~ U, and we must 
set Ci =4/3 to obtain the appropriate Newtonian term 

(7oo = -l + 2C/ + 0(4). 

However, substitution of this value of ci into ((29|) in the same limit gives 

5,, = (l + ;7)(5,y+0(4), 

as previously found by Chiba, Smith and Erickcek [29l |. This result is entirely incompatible with many observations 
of null geodesies (see e.g. [l^) and so we will not consider this case any further. 

Let us now consider case (ii). Here the potentials containing exponentials are expected to be sub-dominant to those 

(2) 

without, on observable length scales, so that the R potential in g^Q is only effective at very small distances. Setting 
ci = 1, and discarding potentials which are exponentially suppressed with regards to others, we then find that the 
metric in this case is given to the appropriate order as 

goo - - 1 + 2C/ - 2C/2 + ^V{pU) - -^V{UR) + 2V(pH) + 4V{pv^) 

9 IStt 

40C2 , , 40 /c7 , , C3 2n 



+ «^(^)-V^^-^yVi^^ + ^'^(^^ (32) 

go^ ^-Iv,- \w, - (33) 



2 2 

(l-f2[/)% (34) 



where '01 and i^ii are defined, as in appendix |B1 by 



_ PV"(X-X0.(X--X") ^3 3 „ 



In the metric above we have retained the term proportional to Zi in the goi components. Although this term contains 
an exponential suppression factor, it is of a different form to Vi and Wi and so cannot necessarily be assumed to be 
negligibly small in comparison to them. 

Another significant difference in the metric above is that the coefficient of the term V{pU) is 64/9, instead of its 
usual value of 4 in GR, and the inclusion of the new potential V{UR). We have already seen that in the limit C2 ^ 
that the V{UR) term reduces to 8TrV{pU)/ci, which is exactly sufficient to recover the GR limit of 4V{pU) in goo- 
The fact that this potential is significant in the limit of small C2, even though it contains an exponential suppression 
term, is due to the factor of C2 in its denominator, as R is present in the integrand. This causes it to approach a 
finite value, instead of zero, when C2 is small. This potential cannot therefore be considered negligible, as it is not 
necessarily exponentially smaller than any other. 

There are three further potentials in goo that are not present in GR: ipi, 1^2 and X{R^). The ipi term contains 
no exponential suppression factor, while the ■02 and X{R^) terms do. However, as before, these potentials are 
included none the less as they are not directly suppressed with respect to any other. The existence of ^"1 in the 
post-Newtonian limit is of particular interest as it is the first new potential that is not exponentially suppressed, with 
obvious significance for constraining the theory with observations. The X{R^) term is also of special interest as it is 
the only term with a dependence on C3. In the limit C2 — ^ this term reduces to 



87rC2 
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which is non-zero when C3 and p 0. This term therefore provides an opportunity to test for deviations from GR at 
the level of in the generating Lagrangian. 

One may now wish to compare the metric obtained to the standard PPN metric, and to read off the relevant 
parameters. The standard PPN metric is given in the present notation by 

-1 + 2U- 2/3C/2 + (27 + 2 + CiWipv^) 

+ 2(37 - 2/3 + 1 + C2)V{pU) + 2(1 + (sWipTl) + 6(7 + C4)V(p) (37) 

-^(47 + 3 + Ci)^.-^(l"Ci)W^. (38) 
(1 + 27C/)(5y (39) 

where /?, 7 and Q are the post-Newtonian parameters, to be set for a particular gravitational theory. We have excluded 
here the preferred location and preferred frame terms, as they are of no relevance for the present study. Comparison 
of the metric ([5^ . ([55)1 and (IM)) with the above allows one to read off the following values 

7-1, Ci=0, C3-O and C4 = 0, 

as in GR. The value of (2 is not so straightforwardly determined, and a naive comparison would yield the result 
(2 = 14/9 instead of the usual value of zero in GR. However, we have seen above that the V{UR) term approaches 
V{pU) for small C2, and gives the GR result in the limit. Care must therefore be taken with the value of this parameter, 
and in the present case it seems more appropriate to consider two contributions towards (^2 - one coming from the 
usual V{pU) term, and the other coming from V{UR). The terms in ([5^ and ([55)1 proportional to tpi, V2, X{R^) and 
Zi are also inadequately accounted for in the PPN metric above. Clearly, new terms are required if these potentials 
are to be included. 



(PPN) 

9oo 



(PPN) _ 

(PPN) ^ 
ilij 



VIII. CONCLUSIONS 



We have determined here the post-Newtonian limit of fourth-order theories of gravity that are analytic functions 
of the Ricci tensor, and that admit Minkowski space as a background. In the Newtonian limit we have recovered the 
well known result that an exponentially suppressed Yukawa potential is present. These deviations from Newton's law 
should be expected to be observed at small distance scales, and a number of experimental efforts have been made to 
find them (see e.g. [30j). These searches have not yet detected any deviations from Newton's law at small distances, 
and so we must consider terms containing exponential factors to be heavily suppressed. 

To determine the post-Newtonian limit of null geodesies we require knowledge of the gij components of the metric 
to an accuracy of 0(2), as is given in equation (p9)) . We again find the well known result that the only correction 
to this term at the post-Newtonian level of accuracy is in the form of a Yukawa potential, which must be heavily 
suppressed. Comparison with the PPN metric ([55]) then gives us that we should expect 7 = 1 for these theories, and 
hence that experiments involving observations of null geodesies should be unable to distinguish them from GR. 

We then proceeded to determine the full post-Newtonian limit to 0(4) in the 500 component, (|?7|) . and to 0(3) in 
the goi component, (|28p . It is found that there exists a large number of new potentials at this order of perturbations. 
Using our knowledge that terms containing exponential factors are heavily suppressed it is possible to neglect a number 
of these potentials, resulting in equations ([32)1 and ((33)) . A comparison of these results with (|37p and p8|) shows that 
the PPN parameters /3, Ci, C3 and (^4 all take the same values in these f{R) theories as they do in GR. Experiments 
which are designed to determine these parameters will therefore be unable to distinguish between the two. However, 
there are differences between ([5^ and ([55]) . and the PPN limit of GR, that may be potentially observable. 

A potentially significant difference with GR is the value of the PPN parameter C2- A direct comparison of ([32]) 
with (|37p appears to yield the result ^2 = 14/9, which is a significant difference from its value of zero in GR. However, 
we know that in the limit that GR is approached the new potential V(UR) in goo makes a contribution that is 
exactly enough to cancel the value of 14/9 above. This strongly suggests that this potential should be included 
when observational constraints are applied. The parameter C2 is usually associated with violations of momentum 
conservation, and the (lack of) self-acceleration of the binary pulsar PSR 1913-1-16 has led to the bound C2 < 4 x 10"'^ 
[3H . Furthermore, it appears likely that observations of the binary system PSR J1738-I-0333 will offer even tighter 
constraints [sl]. However, these previously obtained constraints may not be directly applicable to the current theory. 
Firstly, we know that f{R) theories of gravity covariantly conserve four-momentum exactly (due to their Lagrangian 
formulation). Secondly, these constraints have been imposed in the absence of the V{UR) potential. The extent to 
which these systems are able to offer constraints on f{R) theories when the V{UR) term is included remains to be 
determined. 
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Further opportunity for observationally constraining these theories comes from the other extra potentials in (p2)) 
and ([55)1 : ipi, ip2, X{R'^) and Zi. These potentials are likely to be small, due to suppressing factors of C2/|x — x'p. 
The V'l term is somewhat promising as its suppression is polynomial, and not exponential. The Zi term is also 
interesting as it is has polynomial amplification, as well as exponential suppression. Such a term may be potentially 
observable in experiments that measure vector perturbations, such as gravity probe B [33]. Finally, we will mention 
that the term proportional to X{R'^) allows for the possibility of constraining any R'^ term that may exist in the 
generating Lagrangian. Unlike the R^ term, the R^ term does not first appear in the perturbative expansion as a 
Yukawa potential. In fact, in the limit that the R^ term vanishes the potential X{R^) reduces to p^, which could be 
observable in high density environments. 



APPENDIX A: SOLVING THE t - t EQUATION TO 0(4) 

In order to determine the h^^Q term, it will first be necessary to determine the Ricci scalar to 0(4). This will be 
achieved by solving the trace equation to the appropriate order. In the gauge defined by equations ^ and ([7]), the 
trace of the field equations ^ becomes 

V2i?(4) _ ^i?(4) _ R + ^V^R^ + ^UV'R ~ 2^RV'R - 244i?..,i?.., + &^U,,R,, =—p- 

6C2 2C2 Ci Ci c( cf C2 3C2 

where here we have dropped the superscript on so that the R above should be implicitly assumed to be of 0(2) 
(as can be recognised from the required order of each term) . We now introduce the new potentials 

p'e-/^l"-"'ld3a;' (Al) 



X{Q) ^ I • - d^x', (A2) 

where a prime now denotes a quantity that is a function of x' , so that 

Recognising the relations 

= i ( v^- I VC/ • Vi? + I^Vp • Vf/ + 27rVp • Vi? 
2 \ 6c2 / 3c2 

R^,R^3 - V - — — 1— - + — Vp • Vi? + -J-^R - -^PR 



6c2 / V 2 24c2 / 3c2 144c^ 9c: 
then allows the trace equation above to be integrated, to give 

^/6^ V2c2 cj ' cl 127rc2 V ^4 J Stt 

2X{pU) AXjpR) AX{\7p-VU) AC2X {\I p ■ \I R) X{p) Xjpll) 

3ciC2 3ci 3c^ cf C2 3C2 



(A3) 



We are now sufiiciently equipped to solve the t — t field equation to 0(4), in order to obtain h^QQ ■ Equations ([3]) give 
this as 

_ ^ V^^^) - -V^U' +(lc3- 2^) V'R' - S^V'UR - IS^V^i? 00 
2 ci \2 ci/ ci ci 

+ + C2V'^R'^^^ + 8^[/,y C/,,j - 8^t/,yi?,y - 48^R^,jR^,j 

5 Co n 447r 527rc2 n 

+ -UR - -f - —pU + ---^pR = AnpU + Snpv^ + Uirp 

6 6 3ci 3ci 
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where the trace of ([3]) has been used to ehminate the term proportional to i?^^^ . On recognising 
1. 



TTCl 



^V^UR + -^Vp ■ VU + 2nVp ■ VR - ^pR - 
24c2 3c2 6c2 loci 



pU 



lUcl 



UR 



1ZC2 3C2 



36c^ 



27rci 



pi? 



this can be integrated to 



''00 — 



Cl cf 



16C2 

3cl 



UR — 36— i?,oo 



C2 
— i 



,C2 



V?7| 



8^VU ■ VR 



ISttc 



-V{UR) 



-1 '--1 

-^(^ ) + ^^(pu) - 



48-||Vi?|' 



3c2 ,,,„2n , 64,,, 44c2,,, 40c2, 



47rci ' ' 9cl ^' ' 3cl ^' ' 3cf ^ 



Vf/) 



+ ^V(Vp • VR) + -V(pn) + -V{pv^) + -V{p) + 2^i?(4) 



Cl 



Cl 



Cl 



Cl 



(A4) 



where 



V(Q) 



(A5) 



Equations (|A3[) and (|A4p now specify ftoo to 0(4). 



APPENDIX B: RE-EXPRESSING THE POTENTIALS 

The potentials V(Vp • VU), V(Vp • VR), X{Vp ■ VU) and X{Vp ■ VR) in equation ^ are not in usual PPN form, 
as they are written as functions of gradients of p, U and R. Here we re-express these terms, with the gradients absent, 
as 



V(Vp 




= 47rV(p2) + 


V'l 








V(Vp 


VR) 




-|^V(pi?) + , 


rcT 

V 6c2 


-02 
3C2 


3C2 


A(Vp 


VU) 


= 47rA(p2) 4 










X(Vp 


VR) 


3C2 


- -^X(pi?) + 

6C2 


V'6 
3C2 


f J 


6C2 



18cl 



On substituting these expressions back into ([27| the terms proportional to V(p^) and A(p^) cancel exactly. The new 
potentials, "011 are defined by 



01 = 

02 = 

03 = 

04 = 

05 = 

06 = 



P'p"(x 




x' 


)• 


(x' - 


x") 


|x- 


X 




X' 


-x" 




p'p"(x 




x' 


)• 


(x' - 


x") 


|x- 


X 




X' 


-x" 




pV"(x 




x' 


)• 


(x' - 


x") 




X 




X' 


-x" 




pV"(x 




x' 


)• 


(x'- 


x") 


|x- 


X 




X' 


-x" 




p'p"(x 




x' 


)• 


(x'- 


x") 


|x- 


X 




X' 


-x" 




P'p"(x 




x' 


)• 


(x'- 


x") 


|x- 


X 


p 


X' 


-x" 


3 



x-'|-v/^|x'-x"|^3^,^3^„ 
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_ f p'p"{^-^) ■ (X^ - X") /^|x-x'|- /^|x'-x"| 3 / .3 // 
|X-X'|3|X'-X"|2 ^ "^^''^ 

p'p"{-^ - XQ • (X^ - X'O /^|x-x'|-,/^|x'-x"| 
|X-X'|2|X' -X"|3 
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